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Abstract 

We have analyzed the phenomenon of stochastic resonance in a system driven 
by non Gaussian noises. We have considered both white and colored noises. 
In the latter case we have obtained a consistent Markovian approximation that 
enables us to get quasi-analytical results for the signal-to-noise ratio. As the 
system departs from Gaussian behavior, our main findings are: an enhance- 
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ment of the response together with a marked robustness against noise tuning. 
These remarkable findings are supported by extensive numerical simulations. 
We also discuss the relation with some experiments in sensory systems. 



Stochastic resonance (SR) has attracted considerable interest due, among other aspects, 
to its potential technological applications for optimizing the response to weak external signals 
in nonlinear dynamical systems, as well as to its connection with some biological mechanisms. 
There is a wealth of papers, conference proceedings and reviews on this subject, Ref. [|1| being 
the most recent and complete one, showing the large number of applications in science and 
technology, ranging from paleoclimatology to electronic circuits g], lasers chemical 
systems 0, and the connection with some situations of biological interest (noise- induced 
information flow in sensory neurons in living systems, influence in ion-channel gating, in 
visual perception or in the human blood pressure regulatory system) [^]. Recent works have 
shown a tendency, determined by the possible technological applications, pointing towards 
achieving an enhancement of the system response (by means of the coupling of several SR 
units in what conforms an extended medium P-[TO||), or analyzing the possibility of making 
the system response less dependent on a fine tuning of the noise intensity |JlT|], as well as 
different ways to control the phenomenon [|12|. 

A majority of studies on SR have been made on bistable one-dimensional double-well 
systems, and in almost all cases, with very few exceptions []T3|, the noises are assumed to 
be Gaussian. However, some experimental results in sensory systems, particularly for one 
kind of crayfish |l4| as well as recent results for rat skin |15], offer strong indications that 
the noise source in these systems could be non Gaussian. 

In this letter we analyze the case of SR when the noise source is non Gaussian. We start 
sketching the study of a particular class of Langevin equations having non Gaussian station- 



ary distribution functions [16|, a work based on the generalized thermostatistics proposed 



by Tsallis [17] which has been successfully applied to a wide variety of physical systems 18 
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We consider the following problem 

x = f{x,t)+g{x)7]{t) (1) 

V = ---^VM + -m (2) 
T drj T 

where ^(t) is a Gaussian white noise of zero mean and correlation (^(t)^(t')) = 2 D6(t — t'), 
and Vq{ri) is given by |T6| 

K(^) = ;^(^ln[l + /?(g-l)y], (3) 

where (3 = j^. The function f{x,t) is derived from a potential U{x,t), consisting of a 
double well potential and a linear term modulated by S{t) = F cos{ujt) {f{x,t) = — |^ = 
—Uq + S{t)). For F = this problem corresponds to the case of diffusion in a potential 
Uo{x), induced by rj, a colored non-Gaussian noise. The stationary probability distribution 
for the random variable 7] is given by Pg^irj) = Z^^ 1 + — 1)^ 
normalization factor. Clearly, when g ^ 1 we recover the limit of rj being a Gaussian 
colored noise (Ornstein-Uhlenbeck process). For g < 1 the above distribution cuts-off the 
values of the rj process to | |< [(1 — q)^]^^^"^- It is easy to show that the stationary mean 
value and variance are given, respectively, by {rj) = and (r/^) = [/5(5 — 3g)]~^, the latter 
diverging for all q > 5/3. 

Applying the path-integral formalism to the Langevin equations given by Eqs. (|l]-0), 
and making an adiabatic-like elimination procedure ||19|- pl|] it is possible to arrive at an 



' ^ where Zg is the 



effective Markovian approximation. The specific details are shown elsewhere [^. Such an 
approximation yields the following Fokker-Planck equation (FPE) for the evolution of the 
probability P{x,t) 

dtP{x,t) = -d4A{x)P{x,t)] + ^dl[B{x)P{x,t)]. (4) 

Where 

U' 



and 



B(x) 



D 



rU!^[l + ^{q - + [1 - ^{q - 

The stationary distribution of the FPE in Eq. (^) is thus 



(6) 



B{x) 



exp [— $(x)] 



(7) 



where M is the normahzation factor, and 

Mx) = 2 



A{x') 
B{x') 



dx'. 



(8) 



The indicated FPE and its associated stationary distribution enable us to obtain the mean- 



first-passage-time (MFPT) through a Kramers hke approximation |2^. This quantity is the 
necessary ingredient to work within the two-state model (TST) PJ2^. The MFPT to reach 
X = xq starting from x = a can be obtained from 



T(xo)= r dye'^^y^ T dz B{z)-^e-'^^'\ 

J a J —oo 



(9) 



We focus on polynomial-hke forms for the potential, adopting 
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U{x) = X - y + S{t)x, 



(10) 



setting S{t) = F cos(u;t), and assuming that uj~^ is large compared to the characteristic 
relaxation times in both wells. 

We follow the standard TST approach in order to obtain R, the signal to noise ratio 
(SNR). In our case the result for R, defined as the ratio of the strength of the output signal 
to the broadband noise output evaluated at the signal frequency, is 



R 



F^TT 

2T3 



dT 
dS{t) 



1 2 



(11) 



S(t)=0 



with S{t) the applied signal. 

We present the results for the SNR obtained evaluating Eq. ([IT]). In Fig. la we depict 
the SNR vs. the noise intensity Dc, for a fixed value of the time correlation r and various q, 



while Fig. lb shows the equivalent plot for a fixed value of q and several values of r. However, 
the noise intensity of the non Gaussian process is related to the Gaussian white noise 
D through the scaling Dc = -D[5 — 3g]~^, that was the scaling we used. In the former case 
the general trend is that the maximum of the SNR curve increases when g < 1, i.e. when 
the system departs from the Gaussian behavior. In the latter case the general trend agrees 



with the results for colored Gaussian noises 25 , where it was shown that the increase of 



the correlation time induces a decrease of the maximum of SNR as well as its shift towards 
larger values of the noise intensity. The latter fact is a consequence of the suppression of 
the switching rate with increasing r. Both qualitative trends are confirmed by extensive 
numerical simulations of the system in Eqs. (|l]-0), which has been numerically integrated 
using a stochastic Runge-Kutta-type method |2^. The results were obtained averaging over 
2000 trajectories (5000 trajectories for r = 0). Fig. 2a shows the simulation results for 
the same situation and parameters indicated in Fig. la. Here, in addition to the increase 
of the maximum of the SNR curve for values of g < 1, we also found another remarkable 
aspect that is not reproduced or predicted by our effective Markovian approximation. It is 
the fact that the maximum of the SNR curve flattens for lower values of g, indicating that 
the system, when departing from Gaussian behavior, requires a less flne tuning of the noise 
intensity in order to maximize its response to a weak external signal. Fig. 2b shows the 
simulation results for the same situation and parameters indicated in Fig. lb. Again we 



found an agreement with the general behavior found for colored Gaussian noises 



In order to compare with some experiments where the need to use non-Gaussian noises 



has been suggested, we consider the results obtained for a crayflsh in Ref . [|1J| . A qualitative 
comparison of our results in Figs. 1 and 2 with those given there shows that the agreement 
between theory and experiment does indeed improve by the use of non-Gaussian noises 
corresponding to g < 1. However, it is known that rather than applying our simple one- 
variable model, to the crayflsh neural system under consideration it is better to use the 



celebrated FitzHugh-Nagumo one. In the same way as in Ref. |T^, we deflne the model 
system through the set of equations 



Tyi) = v{v — 0.5)(1 — v) — 

TwW = V — w + ecos(fit) + ri{t), (12) 

where v{t) is the variable associated to the action potential (in an activator-inhibitor model, 
it corresponds to the activator variable), and w{t) is the recovery (inhibitor) variable. The 
characteristic times for such variables are, respectively, and Tw The recovery variable is 
submitted to a periodic signal (whose (small) amplitude we denote by e and its frequency by 
Q) and a noise source ri{t). This random variable, now chosen as white (that is with r = 0), 
is distributed accordingly to the stationary pdf of Eq. (13). We have integrated numerically 



Eqs. (|T2[) using a stochastic Runge-Kutta-type method with a time step of At = 1ms. 

The time series v{t) was converted into a spike train (mimicking what happens in the 
nervous system). One spike occurs when v{t) exceeds the threshold potential be- Each spike 
was modeled as a square wave of height Vs = IV, and duration tg = 3ms. It is not possible 
to have two successive spikes with a lag smaller than tg- The SNR was evaluated following 
the usual procedure on this modified signal. 

Along all the simulations, we fixed = 10~^s and = lO^^s. While we considered the 
signal as having frequency Q = 55Hz and amplitude e = 0.03V. All these values agree with 
those from Ref. |l^, however their results are a particular case of ours, i.e. q = 1, for the 



noise source. 

We have compared our simulations for R, the SNR, as a function of noise standard 
deviation for different values of q and be, the threshold potential, with the experimental data 
1^ . Here again we have scaled the noise intensity according to Dn = D[5 — 3g]^^, where 



m 



Dn is the intensity of the Non Gaussian noise while /^corresponds to the pdf parameter. 
Figure 3 shows the results obtained for be = 0.18V. We see that the (Gaussian) case q = I 
underestimates the measured SNR values for the whole range of noise intensities, while 
q = 0.75 exhibits an overestimation for large noise values. However, if we adopt q = 0.47 
the system exhibits a nice agreement with experimental data both near the SNR peak and 
for large values of noise. When the threshold value is fixed at be = 0.15V as in (results 
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not shown here), the case q = 1 corresponds exactly to the theoretical fit in [0. From 



the comparison with g 7^ 1 it becomes apparent that, although the fit is acceptable near 
the R maximum, for large noise intensities the q = 1 case underestimate the values of R. 
By setting q = 0.75, we find a much better fit even in the region of large noises, and when 
q = 0.30 and for large values of the noise standard deviation, the results exhibit values above 
the measured ones. 

It it worth remarking here that the experimental behavior in the low noise limit is not 
reproduced by any model. This is due to the spontaneous firing of neurons in the real 



crayfish neural system 14 



we 



Summarizing, motivated by some experimental results in sensory systems |ll^ , p!5 
have analyzed the problem of SR when the noise source is non Gaussian. We have chosen 
a non Gaussian noise source (white or colored) with a probability distribution based on the 
generalized thermostatistics fl^ . In the colored case and making use of a path integral ap- 
proach, we have obtained an effective Markovian approximation that allows us to get some 
analytical results. In addition, we have performed exhaustive numerical simulations. Even 
though the agreement between theory and numerical simulations is only partial and qual- 
itative, the effective Markovian approximation turns out to be extremely useful to predict 
(qualitatively) general trends in the behavior of the system under study. 

Our numerical and theoretical results indicate that: (i) for a fixed value of r, the max- 
imum value of the SNR increases with decreasing q; (ii) for a given value of q, the optimal 
noise intensity (that one maximizing SNR) decreases with q and its value is approximately 
independent of r; (iii) for a fixed value of the noise intensity, the optimal value of q is 
independent of r and in general it turns out that q ^ 1- 

As we depart from Gaussian behavior (with g < 1), the SNR shows two main aspects: 
firstly its maximum as a function of the noise intensity increases, secondly it becomes less 
dependent on the precise value of the noise intensity. Both aspects are of great relevance for 
technological applications However, as was indicated in Ref. [|r^, non Gaussian noises 
could be an intrinsic characteristic in biological systems, particularly in sensory systems 



||7| JT^JT5[] . In addition to the increase in the response (SNR), the reduction in the need of 
tuning a precise value of the noise intensity is of particular relevance both in technology and 
in order to understand how a biological system can exploit this phenomenon. As an example 
of this, simulations of a FHN model with non Gaussian noise shows a behavior in better 
agreement with the experimental results in a crayfish sensory system than simulations made 
with a Gaussian noise source. The present results indicate that the noise model used here 
offers an adequate framework to analyze such a problem. 

Acknowledgments: The authors thank C. Tsallis for useful discussions and V. Grunfeld 
for a revision of the manuscript. HSW acknowledges support from CONICET, Argentine 
agency. RT acknowledges financial support from DGES, project numbers BFM2000-1108 
and PB97-0141-C02-01. 



8 



REFERENCES 



[1] L. Gammaitoni, P. Hanggi, P. Jung and F. Marchesoni, Rev. Mod. Phys. 70 223 (1998). 

[2] R. Benzi, A. Sutera and A. Vulpiani, J. Phys. A 14, L453 (1981); R. Benzi, G. Parisi, 
A. Sutera and A. Vulpiani, Tellus 34, 10 (1982). 

[3] C. Nicolis, Tellus 34, 1 (1982). 

[4] S. Fauve and F. Heslot, Phys. Lett. A 97, 5 (1983); R.N. Mantegna and B. Spagnolo, 
Phys. Rev. E 49, R1792 (1994); V.S. Anishchenko,M.A. Safonova and L.O. Chua, Int. 
J. Bifurcation and Chaos Apph. Sci. Eng. 4, 441 (1994). 

[5] J.M. lanneUi, A. Yariv, T.R. Chen and Y.H. Zhuang, Appl. Phys. Lett. 65, 1983 (1994); 
A. Simon and A. Libchaber, Phys. Rev. Lett. 69, 3375 (1992). 

[6] A. Guderian, G. Dechert, K. Zeyer and F. Schneider; J. Phys. Chem. 100, 4437 (1996); 
V. Petrov, Q.Ouyang and H.L. Swinney Nature 388, 655 (1997). 

[7] J. K. Douglas, L. Wilkens, E. Pantazelou and F. Moss, Nature 365, 337 (1993); S. M. 
Bezrukov and I. Vodyanoy, Nature 378, 362 (1995); E. Simonotto, M. Riani, C. Seife, 
M. Roberts, J. Twitty and F. Moss, Phys. Rev. Lett. 78, 1186 (1997); I. Hidaka, D. 
Nozaki and Y. Yamamoto, Phys. Rev. Lett 85, 3740 (2000). 

[8] J.F. Lindner, B.K. Meadows, W.L. Ditto, M.E. Inchiosa and A. Bulsara, Phys. Rev. 
Lett. 75, 3 (1995), and Phys. Rev. E 53, 2081 (1996). 

[9] F. Castelpoggi and H.S. Wio, Europhysics Letters 38, 91 (1997); F. Castelpoggi and H. 
S. Wio, Phys. Rev. E 57, 5112 (1998); M. Kuperman, H.S. Wio, G. Izus, R. Deza and 
F. Castelpoggi, Physica A 257, 275 (1998). 

[10] S. Bouzat and H. S. Wio, Phys. Rev. E 59, 5142 (1999). 

[11] C.J. Tessone and H.S. Wio, Mod.Phys.Lett.B 12, 1195-1202 (1998); C.J. Tessone, H.S. 
Wio and P.Hanggi, Physical Review E 62, 4623 (2000). 

9 



[12] L. Gammaitoni, M. Locher, A. Bulsara, P. Hanggi, J. Neff, K. Wiesenfeld, W. Ditto 
and M. E. Inchiosa, Phys. Rev. Lett. 82, 4574 (1999). 

[13] D. Nozaki, J. Collins and Y. Yamamoto, Phys. Rev. E 60, 4637 (1999). 

[14] K. Wiesenfeld, D. Pierson, E. Pantazelou, Ch. Dames and F. Moss, Phys. Rev. Lett. 
72, 2125 (1994). 

[15] D. Nozaki, D.J. Mar, P.Grigg and J.J.Collins, Phys. Rev. Lett. 72, 2125 (1999). 

[16] L. Borland, Phys. Lett. A 245, 67 (1998); L. Borland Phys. Rev. E 57, 6634 (1998). 

[17] E. M. F. Curado and C. Tsallis, J. Phys. A 24, L69 (1991); ibidem 24, 3187 (1991), 
ibidem 25, 1019 (1992). 

[18] A.R. Plastino and A. Plastino, Phys. Lett. A 174, 384 (1993); D.H. Zanette and P.A. 
Alemany, Phys. Rev. Lett. 75, 366 (1995); B.M. Boghosian, Phys. Rev. E 53, 4754 
(1996); C. Tsallis and A.M.C. de Souza, Phys. Lett. A 235, 444 (1997), G. Drazer, 
H.S.Wio and C. Tsallis, Phys. Rev. E 61, 1417 (2000); R. Salazar and R. Toral, Physica 
A290, 159 (2001). 

[19] P. Colet, H.S. Wio and M. San Miguel, Phys Rev A39, 6094 (1989); H.S.Wio, P.Colet, 
L.Pesquera, M.A.Rodriguez and M.San Miguel, Phys.Rev. A 40, 7312 (1989). 

[20] F.Castro, H.S.Wio and G.Abramson, Phys.Rev.E 52, 159 (1995); F. Castro, A. Sanchez 
and H. S. Wio, Phys. Rev. Lett. 75, 1691 (1995). 

[21] H.S. Wio, Application of Path Integration to Stochastic Processes: An Introduction, in 
Fundamentals and Applications of Complex Systems, Ed. G.Zgrablich (Nueva Ed. Univ., 
U.N.San Luis, 1999), pg.253. 

[22] M. A. Fuentes, Raul Toral and Horacio S. Wio, Effective Markovian approximation for 
Non Gaussian Noises: a path integral approach, submitted to Physica A. 

[23] P. Hanggi, P. Talkner and M. Borkovec, Rev. Mod. Phys, 62, 251 (1990). 



10 



[24] B. McNamara and K. Wiesenfeld, Phys. Rev. A 39, 4854 (1989). 

[25] L. Gammaitoni, E. Manichella-Saetta, S. Santucci, F. Marchesoni and C. Presilla, Phys. 
Rev. A 40, 2114 (1989). 

[26] P. Jung and P. Hanggi, Europhys. Lett. 8, 505 (1989); P. Jung and P. Hanggi, Phys. 
Rev. A 44 8032 (1991). 

[27] M. San Miguel and R. Toral, Stochastic Effects in Physical Systems, in Instabilities and 
Nonequilibrium Structures VI, E. Tirapegui, J. Martinez and R. Tiemann, eds., (Kluwer 
Ac. Pub., 2000). 



Figure 1: Theoretical value of the SNR, R, vs. the noise intensity D, for: (a) r = 0.1 and 
the following values of g = 0.25, 0.75, 1.0, 1.25 (continuous line, from top to bottom), 
(b) q — 0.75 and the following values of r = 0.25,0.75, 1.5 (broken hne from top to 
bottom). We adopted F = 0.1 and uj = 0.1. 

Figure 2: Simulation results of the SNR, R, vs. the noise intensity D, for (a) r = 0.1 and 
the following values of g = 0.25, 0.75, 1.0, 1.25 (from top to bottom), (b) q = 0.75 and 
the following values of r = 0.25, 0.75, 1.5 (from top to bottom). We adopted F — 0.1 
and u! — 0.1. 



Figure 3: The SNR, R, for the FHN model is plotted as a function of V^D (the noise 
standard deviation). The different curves show the results for different g— values: 
g = 1 (•), g = 0.75 (0), g = 0.47 (V). In all the cases we fixed be = 0.18V, simulation 
results are compared with the experimental data (□). 
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